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Effects of randomness on interacting fermionic systems 
in one dimension are investigated by quantum Monte-Carlo 
techniques. At first, interacting spinless fermions are studied 
whose ground state shows charge ordering. Quantum phase 
transition due to randomness is observed associated with the 
collapse of the charge ordering. We also treat random Hub- 
bard model focusing on the Mott gap. Although the random- 
ness closes the Mott gap and low- lying states are created, 
which is observed in the charge compressibility, no (quasi-) 
Fermi surface singularity is formed. It implies localized na- 
ture of the low-lying states. 
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Effects of Coulomb interaction and randomness in 
fermionic systems have been studied separately for a long 
time. Their combined effects, however, remain unsolved 
especially in the strong-coupling regime. The strong 
Coulomb interaction can cause the Mott transition and 
lead the metallic system to the Mott insulator. Also the 
randomness can make the electrons localized, which is 
the Anderson localization. Although both are metal- 
insulator transitions, the physical characters are quite 
different. For example, in the Mott insulator, the ground 
state has a finite charge gap due to the Coulomb inter- 
action. On the other hand, the insulator due to the An- 
derson localization does not necessarily have a charge 
gap. Therefore the competition between the Coulomb 
interaction and the randomness is interesting from the- 
oretical point of view 0-f§ and it may shed new light 
on some real materials. In this paper, we systematically 
investigate one-dimensional fermionic systems with both 
Coulomb interaction and randomness by quantum Monte 
Carlo (QMC) techniques |-(|. Although the conven- 
tional wisdom is that there is no quantum phase tran- 
sition driven by randomness in one and two dimensions 
we confirmed the existence of quantum phase tran- 
sitions due to combined effect of Coulomb interaction 
and randomness. In this paper, at first, we focus on 
the charge density wave (CDW) ordering in interacting 
spinless fermions with randomness. Next, in order to 
understand low-lying excitations of random Mott insula- 
tor, we study the charge compressibility of random Hub- 
bard model in the half-filled sector. Also the one-particle 
Green's function is investigated to discuss nature of the 
low-lying excitations. 

Random Spinless Fermions- To begin with, we focus 
on CDW structure factor Ocdw of interacting spinless 
fermions with randomness in one dimension. The Hamil- 
tonian is given by 



Ti-sf = 7~(v + T~i w , (1) 

H v = -tJ2 (cUj + c] Ci ) + VJ2 n inj , (2) 
(id) (id) 

7-L w =22 w i n h ( 3 ) 
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where c\(ci) creates(annihilates) a spinless fermion at the 
i-th site, ni = ejej and (i,j) is a nearest- neighbor link. 
TL W denotes site randomness where Wi's take W or — W 
at random. We set t = 1 as an energy unit and use open 
boundary condition. In the following, the system is set 
to be half-filled i.e. N e = L/2, where N e is number of 
the electrons and L is the system size. The pure system 
without randomness is one of the typical examples which 
show quantum critical phenomena. When the interaction 
is weak ( 0<V/t<2 ), the system is metallic and it be- 
longs to a universality class of the Tomonaga-Luttinger 
liquid. A metal-insulator transition accompanied with 
a charge ordering occurs at V/t = 2. For strong inter- 
action (V/t > 2), the charge degree of freedom becomes 
frozen and the energy gap opens. We focus on the charge- 
ordered phase {V/t > 2) where the charge gap is finite 
due to the nearest-neighbor interaction, and study the 
effect of randomness. 

In order to study the charge ordering when the ran- 
domness is included, we calculate the CDW structure 
factor with momentum 7r defined by 

Ocdw(T) = ^ E - 1/2)K - 1/2)), (4) 

where B is set to be a subset of the lattice to avoid a 
boundary effect, B = [L/4, 3L/4] and j}£> is a number 
of sites in B. Without randomness, Ocdw is strongly 
enhanced and diverging toward zero temperature which 
corresponds to charge ordering in the ground state due to 
Coulomb interaction (V/t > 2). On the other hand, when 
sufficiently strong randomness is included, the charge or- 
dering is expected to collapse and the diverging behavior 
of Ocdw to vanish. In order to obtain the Ocdw for 
each realization of randomness, world-line QMC method 
is employed ||. In Fig.l, temperature dependence of 
the Ocdw is shown for several strength of randomness 
where simulations are performed in the half-filled sec- 
tor with typical system size L = 128 and V/t = 3. 
Without randomness, since the charge ordering occurs 
at zero temperature, the structure factor shows diverg- 
ing behavior toward zero temperature. For weak ran- 
domness, the structure factor still has similar behavior 
down to the temperature we studied. It means (quasi-) 
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long-range order in the ground state or long localization 
length beyond the available system size. On the other 
hand, when sufficiently strong randomness is included, 
the temperature dependence of the Ccdw shows qual- 
itatively different behavior. The rapid enhancement at 
low temperature vanishes. It implies that the charge or- 
dering completely fades out due to randomness. In order 
to investigate detailed nature of the transition, we study 
how the structure factor depends on the randomness at 
T = Q.2t. As shown in the inset of Fig.l, sudden decrease 
of Ocdw at finite strength of randomness is observed. It 
implies that the collapse of the charge gap occurs at fi- 
nite strength of randomness, which is the order of the 
charge gap, although we can not exclude a possibility 
that the transition is the Kosterlitz-Thouless (KT) type 
at infinitesimally weak randomness. 

Random Hubbard model- Next, in order to study effect 
of randomness on the Mott insulator, let us consider half- 
filled sector of one-dimensional random Hubbard model. 
The Hamiltonian is given by 



n u = -* £ (4 



(5) 
(6) 
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where t is the nearest-neighbor hopping amplitude and U 
is the on-site Coulomb interaction. TL W denotes random 
potentials and Wi 's are taken at random from the interval 
[— w,w]. We treat the system in a grand canonical en- 
semble with the chemical potential \i. The boundary con- 
dition is periodic. In the absence of randomness, an in- 
finitesimal interaction U causes a charge gap (Mott gap) 
Eg at half filling (fj, = U/2). The charge gap E g is ex- 
ponentially small in the weak-coupling region (U/t <C 1) 
and linear in U in the strong-coupling region (U It ^> 1), 
for example, estimated E g ~ 1.3t for U/t = 4 M. Here 
we shall discuss effect of randomness on the Mott insu- 
lator (see also Ref. 
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To obtain approximation-free results, we use a finite- 
temperature auxiliary- field QMC method HQ]. Since we 
use the grand canonical ensemble, there is a finite charge 
fluctuation which is crucial for the knowledge of low- lying 
excitations. Although the random potential breaks the 
particle-hole symmetry for each realization of random- 
ness, half-filling condition is recovered after averaging 
over different realizations of randomness. Our simula- 
tions are performed with system size L = 36 at U/t = 4. 
Severe finite-size effect due to energy discretization is ob- 
served in low-tempcrature region (lower than T ~ 0.2t 
for L = 36) and the data for that region are not shown. 
Moreover, since the particle- hole symmetry is broken for 
each realization of randomness, negative-sign problem oc- 
curs in general. However, in the parameter region we 
investigated, it is not serious and the data are obtained 
with sufficient accuracy. To investigate the change of low- 



lying excitations due to randomness, we calculate charge 
compressibility n defined by 
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It measures fluctuation in the charge sector and shows 
thermally-activated behavior when the system has a fi- 
nite charge gap. On the other hand, in the absence of 
the charge gap, the charge compressibility is expected to 
be finite due to the low-lying excitations. For example, 
in the non-interacting case, k(T — 0) is equal to the den- 
sity of states at the Fermi energy. Fig. 2 shows temper- 
ature dependence of the charge compressibility n for dif- 
ferent strength of randomness. Although a snap-shot for 
a particular realization of randomness is shown in Fig. 2, 
the simulations for other realizations were also performed 
and we confirmed that the global feature does not depend 
on each realization. Without randomness, k decreases 
toward zero and shows thermally-activated behavior as 
the temperature is lowered indicating the existence of a 
charge gap. This is a typical feature of the Mott insu- 
lator. As randomness is turned on, the enhancement in 
k is observed for all temperature. For weak randomness, 
w/t<1.5, although k is enhanced, it still seems to show 
thermally-activated behavior as in the pure case. On the 
other hand, with strong randomness w/t>l.b, there is no 
tendency for k to decrease down to the lowest tempera- 
ture we studied and it seems that n is finite at T = 0. It 
suggests the collapse of the charge gap. These behaviors 
imply that sufficiently strong randomness, which is the 
order of the charge gap, takes the system away from the 
Mott insulator |5]4|]. This is a quantum phase transition 
of the Mott insulator driven by randomness. Although 
it has been confirmed that low-lying excitations are cre- 
ated in the Mott gap with sufficiently strong random- 
ness, the nature of the low-lying excitations is crucial 
to understand the phase transition. In order to study 
the low-lying excitations, we calculate the one-particle 
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Green's function Gij a — 
in metallic states is clear in the momentum-space repre- 
sentation, which is the momentum distribution function 
n a (k). Metallic nature of the system is reflected by the 
singularity of n a (k) at k = kp. For example, it is a step 
function in the Fermi liquid (Fermi surface) and, even in 
the Tomonaga-Luttinger liquid (doped Hubbard chain), 
the derivative of n a {k) is diverging at k — kp (quasi- 
Fermi surface). On the other hand, in the Mott Insulator 
(half- filled Hubbard chain), the derivative is finite. We 
may summarize that the formation of (quasi-) Fermi sur- 
face, as temperature is lowered, implies metallic nature 
of the system. Fig. 3 shows momentum distribution func- 
tion n a (k) at T = Q.2t for several randomness strength. 
By performing calculations with different temperatures, 
we confirmed that the temperature is sufficiently low that 
the results reflect the ground-state properties. Although 
the charge gap is closed for w/t>l .5 as discussed above, 
the momentum distribution function is almost unchanged 
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and no singularity appears at Uf even for w/t — 2.0. It 
is in contrast to the transition from the Mott insulator 
to the Luttinger liquid in the pure Hubbard model when 
the system is doped. It suggests that the low-lying exci- 
tations created by randomness are localized. 

In summary, we have studied effects of randomness in 
the interacting fermionic system with charge gap (Mott 
gap). The strong randomness closes the charge gap and 
low-lying states are created, while (quasi-) Fermi surface 
is not formed. It implies that the transition is an insula- 
tor to an insulator transition. The transition point may 
be continuously connected to the metal-insulator transi- 
tion point in the pure system, then it gives a new example 
of quantum phase transitions in low-dimensional random 
fermionic systems. These phenomena may give a clue for 
the understanding of several aspects in quantum phase 
transitions in low-dimensional random systems. 
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FIG. 1. The CDW structure factor as a function of tem- 
perature for a particular realization of randomness (L — 128 
and V/t = 3). For weak randomness (W/t = 0.1), the struc- 
ture factor shows diverging behavior down to the tempera- 
ture we studied. On the other hand, for strong randomness 
(W/t = 3), the divergence vanishes. The inset shows how the 
structure factor at T/t = 0.2 depends on randomness, where 
L — 64, V/t — 3 and average over 30 — 40 realizations of 
randomness is performed. The line is guide for eyes. 

FIG. 2. Temperature dependence of the charge com- 
pressibility k for a particular realization of randomness 
(L = 36 and U/t — 4). Without randomness, k shows ther- 
mally-activated behavior and decreases toward zero at low 
temperature indicating the existence of charge gap. With 
weak randomness, w < w c (w c /t^l.5), k still shows ther- 
mally-activated behavior. On the other hand, for strong ran- 
domness, w > w c , « does not decrease down to the tempera- 
ture we studied. 

FIG. 3. The momentum distribution function (L — 36 and 
U/t = 4). Average over about 10 realizations of randomness 
is performed. There is no singularity at k — &f even for 
w/t = 2.0 where the charge gap is closed as shown in the 
charge compressibility 
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